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Optical variability of the ultraluminous X-ray source NGC 1313 X-2 
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We analyzed the longest phase-connected photometric dataset available for NGC 1313 X-2, looking for the ~6 day mod- 
ulation reported by Liu et al. J2009I I. The folded B band light curve shows a 6 days periodicity with a significance slightly 
larger than Scr. The low statistical significance of this modulation, along with the lack of detection in the V band, make its 
identification uncertain. 
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1 Introduction 

NGC 1313 X-2 is located in the outskirts of the barred spi- 
ral galaxy NGC 1313 at a distance of 3.7-4.27 Mpc (TuUy 
[19881 Mendez et al. l2002l Rizzi et al. lIOOTll . Its observed 
X-ray luminosity varies between a few x 10^^ erg/s and ~ 
10^° erg/s in the 0.3-10 keV band (Feng & Kaaret|2006ll. 
The source has been extensively studied in the X-ray and 
optical bands (e.g. Mucciarelli et al. 120071 Grise et al. l2008l) . 
It belongs to a handful of ultraluminous X-ray sources (ULXs) 
clearly associated to stellar optical counterparts (e.g. Liu et 
al. 2004; Kaaret et aL |2004l Mucciarelli et al. |2055l Soria et 
al. 120051) . These sources appear almost ubiquitously hosted 
in young stellar environments (e.g. Pakull et al. 2006, Ram- 
sey et al 120061 Liu et al. I20071 l and have properties con- 
sistent with those of young, massive stars. However, some 
ULXs appear to be associated to older stellar populations 
and one possible later type stellar counterpart is now known 

al- 



Table 1 Log of the VLTh-FORS 1 and HSTH-WFPC2 pho- 
tometric observations of NGC 1313 X-2 



(Feng & Kaaret '2008; Roberts, Levan & Goad '200811 
though its spectral classification may be affected by sig- 
nificant galactic and extra-galactic reddening (Grise et al. 
|2006|. In the case of NGC 1313 X-2, a single optical coun- 
terpart has been identified through a chain of efforts (Zampieri 
et al. 120041 Mucciarelli et al . 120051120071 PakuII et aL l2056l 
Liu et al.'2007; Grise et al. 1 20081) . This star has an extinction- 
corrected absolute magnitude Mb ~ —4.5 mag and colors 
{B ~ y)o - -0.15 mag and {V -I)o ^ -0.16 mag (Muc- 
ciarelli et al. 120071 Grise et al. 120081 ). consistent with a B 
spectral type. 

Liu et al. ( 120091 ) found a possible periodicity of 6.12± 
0.16 days in the B band light curve of the optical counter- 
part of NGC 1313 X-2, that was interpreted as the orbital 
period of the binary system. Three cycles were detected in 



Obs. 


Date 


MJD 


Exposure 
(s) 


Instr. 


1 


2007 


-10-21 


54394.362587 


242x2 


FORSl 


2 


2007 


-11-15 


54419.220472 


242x2 


FORSl 


3 


2007 


-11-15 


54419.277624 


242x2 


FORSl 


4 


2007 


-11-16 


54420.255469 


242x2 


FORSl 


5 


2007 


-12-06 


54440.075805 


242x2 


FORSl 


6 


2007 


-12-10 


54444.173203 


242x2 


FORSl 


7 


2007 


-12-14 


54448.121299 


242x2 


FORSl 


8 


2008 


-01-31 


54496.073081 


242x2 


FORSl 


9 


2008 


-03-02 


54527.059338 


242x2 


FORSl 


10 


2008 


03-05 


54530.053776 


242x2 


FORSl 


11 


2008 


03-08 


54533.056052 


242x2 


FORSl 


12 


2008 


05-21 


54607.911122 


500x2 


WFPC2 


13 


2008 


05-22 


54608.043761 


500x2 


WFPC2 


14 


2008 


05-23 


54609.042372 


500x2 


WFPC2 


15 


2008 


05-24 


54610.040983 


500x2 


WFPC2 


16 


2008 


-05-25 


54611.038900 


500x2 


WFPC2 


17 


2008 


05-26 


54612.104178 


500x2 


WFPC2 


18 


2008 


05-27 


54613.179178 


500x2 


WFPC2 


19 


2008 


05-28 


54614.178483 


500x2 


WFPC2 


20 


2008 


05-29 


54615.177789 


500x2 


WFPC2 


21 


2008 


-05-30 


54616.097233 


500x2 


WFPC2 


22 


2008 


-05-31 


54617.109039 


500x2 


WFPC2 


23 


2008 


06-01 


54618.692372 


500x2 


WFPC2 


24 


2008 


06-02 


54619.173622 


500x2 


WFPC2 


25 


2008 


06-03 


54620.105567 


500x2 


WFPC2 


26 


2008 


06-04 


54621.171539 


500x2 


WFPC2 


27 


2008 


-06-05 


54622.104178 


500x2 


WFPC2 


28 


2008 


06-06 


54623.102789 


500x2 


WFPC2 


29 


2008 


06-07 


54624.034733 


500x2 


WFPC2 


30 


2008 


06-08 


54625.100706 


500x2 


WFPC2 


31 


2008 


06-09 


54626.611817 


500x2 


WFPC2 
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Fig. 1 Joint VLT+HST light curve of NGC 1313 X-2 in the B band. The squares are the VLT+FORSl data, while the 
circles represent the HST+WFPC2 observations. The magnitudes are the difference with those of a reference field star The 
data cover a period of ^ 7.5 months. The solid (cyan) line is the best fitting sinusoid with a period P = 6.0 days. 



the B band, while no modulation was found in V. Previ- 
ous studies carried out on the available HST and VLT ob- 
servations led to negative results (Grise et al. '2008). More 
recently, lack of significant photometric variability on a new 
sequence of VLT observations has been reported by Grise et 
al. (l2009l l. 

Here we present a preliminary reanalysis of the joint 
VLTh-FORSI and HSTH-WFPC2 photometric observations 
of NGC 1313 X-2 obtained during the years 2007-2008, 
with the aim of clarifying the statistical significance of the 
orbital periodicity identified by Liu et al. ( I2009I I. Further de- 
tails will be reported in a separate paper (Zampieri et al., in 
preparation). 

2 VLT and HST observations 

NGC 1313 X-2 was observed with VLTh-FORSI between 
October 2007 and March 2008 (11 epochs; Grise et al. |2009] l 
and with HSTH-WFPC2 between May and June 2008 (20 
epochs; Liu et al. l2009l l. The quality of the HST images is 
fair, despite the degradation of the central PC chip. A log of 



the observations is reported in Table 1 . We re-analyzed the 
whole dataset in a homogeneous way, looking for the ^6 
day periodicity reported by Liu et al. (|2009l l. 

After performing standard image reduction in the IRAF 
environment, the two exposures taken a few minutes apart 
each night were combined together and cleaned for cos- 
mic rays. To accurately photometer the objects, we used 
AIDA (Astronomical Image Decomposition and Analysis; 
Uslenghi & Falomo l2008b . an IDL-based package originally 
designed to perform two-dimensional PSF model fitting of 
quasar images. For the analysis of the WFPC2 exposures, 
we loaded into AIDA the appropriate PSF simulated with 
Tiny Tim v. 6.fl 

Analyzing VLT and HST measurements together requires 
attention to be paid to the systematic differences between 
the two photometric systems. We then first converted the 
HST instrumental magnitudes to the standard U BVRI pho- 
tometric system using the updated transformation equations 
and coefficients published in Dolphin ( 120091 1 for the appro- 



http://www.stsci.edu/software/tinytim/tinytim.html 
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Fig. 2 Binned light curve (6 bins) of the B band VLT+HST dataset of NGC 1 3 1 3 X-2, folded over the best estimate of 
the period. The best fitting sinusoid is also shown. 



priate instrumental gain (which is equal to 7 in our casefl 
The color correction term was computed adopting the {B — 
V) color reported in Mucciarelli et al. (2005 ). In spite of 
this, residual systematic differences between the two pho- 
tometric systems might still be present and affect our mea- 
surements. In particular the color correction term is sensi- 
tive to the overall bandpass (telescope plus atmospheric re- 
sponse) of the instruments. Further attention should be paid 
to minimize the effects of possible absolute calibration un- 
certainties. We then decided to perform differential photom- 
etry of the target with respect to a nearby field star (star D 
in Zampieri et al. 120041 ). located on the same chip in both 
instruments. The reference star is brighter than the target 
and has a low root mean square variability (~ 0.05 mag 
in the VLT and ^ 0.02 mag in the HST exposures). For 
similar reasons, during all the observations performed with 
HSTH-WFPC2, the field was always oriented in the same di- 
rection and the target and reference star were always located 
on the same position on the central PC chip. 

Figure [T] shows the B band light curve of NGC 1313 X- 
2 obtained in this way. The data show clear short term (^ 1 
day) variability, likely due to X-ray irradiation. As can be 
seen from Figure [1] the VLT data have much smaller root 
mean square variability (^0.04 mag) than the HST ones 
(~0.09 mag). Thus, there are different levels of optical ac- 
tivity whose origin is unclear and is under investigation at 
present (Zampieri et al., in preparation). Superimposed on 
the short term stochastic variability, the HST dataset shows 
also an approximately sinusoidal modulation with a period 
of 6 days (Figure[T]|. 



3 Results 



Following Liu et al. (12009'), we fitted all the VLTh-HST 
datasets (31 epochs) with a sinusoid: 

AB ^ B + Asm{2n{t~ti)/P + <j)), (1) 

where A, P and ip are the amplitude, period and phase, 
respectively, ti = 54390 is a reference epoch and B = 
< AB >—< Bd — Btarget >= —3.198 is the average 
differential magnitude. The best fitting parameters are: P = 



6.01 



+0.17 
-0.01 



days, A 



0.09 



+0.02 
-0.02 



mag. 



= 94' 



+ 13" 



The 



http ://purcell . as . arizona.edu/wfpc2 _ calib/ 



value of the amplitude and period are in agreement with 
those reported by Liu et al. ( 120091 ). The error on the pe- 
riod is such that AP{T/P) <, 6 days, where T - 200 
days is the interval between the first and last observation 
(sampling interval). This is not larger than P, indicating that 
the two datasets can be phase connected. Although the VLT 
data alone do not show evidence of periodicity, they appear 
to be consistent with the sinusoidal modulation observed in 
the HST observations. A comparison of the VLT data alone 
with the best fit from the HST data shows that the former are 
statistically consistent with the HST fit (reduced chi-square 
Xr = 0.42). 

We reanalyzed also the V band HST+WFPC2 images 
and found no significant periodic variability. If the observed 
modulation is caused mostly by X-ray irradiation, the am- 
plitude in the V band is expected to be smaller than that in 
the B band (Patruno & Zampieri 120 101 ). A sinusoidal mod- 
ulation with the same period and phase obtained from the 
fit of the B band data and an amplitude 0.04 mag is con- 
sistent, within the errors, with the data (reduced chi-square 

^ 1.3). 

The statistical significance of the B band modulation 
was tested performing a Lomb-Scargle periodogram anal- 
ysis of all the observations. We found that the modulation 
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is significant only at the ^ l.dcr level. Binning the light 
curve in 6 bin intervals and performing an epoch folding 
period search, the 6 days modulation is recovered with a 
significance slightly larger than Su (see Figure|2]l. Although 
the binned light curve suggests that the periodicity may be 
there, the low statistical significance of the B band modula- 
tion, along with the lack of detection in the V band, make 
its identification uncertain. A dedicated photometric moni- 
toring campaign under homogeneous observing conditions 
to minimize systematic uncertainties are needed to confirm 
it. 
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